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Sylvester [20] asks for the circle of smallest radius enclosing a given set of finite points in the 



^vq | plane. An extension of the smallest enclosing circle problem called the smallest intersecting ball 



problem was considered in [17, 18]: given a finite number of nonempty closed subsets of a normed 
space, find a ball with the smallest radius that intersects all of the sets. In this paper we initiate the 
study of minimal time functions generated by unbounded dynamics and discuss their applications to 
extensions of the smallest intersecting ball problem. This approach continues our effort in applying 
convex and nonsmooth analysis to the well-established field of facility location. 
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1 Introduction and Preliminaries 

Given a finite number of points on the plane, the sn 

the smallest circle that encloses all of the points. This celebrated problem was introduced in 



(N 

J> ■ Given a finite number of points on the plane, the smallest enclosing circle problem asks for 

OO 
00 . 

\Q . the 19th century by the English mathematician J. J. Sylvester [20] and is a simple example of 

continuous facility location problems. The reader is referred to [7, 19, 21] and the references 
therein for more recent results involving this problem in both theoretical and numerical 
aspects. 

Let X be a normed space and let F C X be a nonempty closed convex subset of X. 
Given nonempty closed target sets Qj for i = 1 , . . . , m and a nonempty closed constraint 
set (!CI, our goal in this paper is to study the following extended version of the smallest 
intersecting ball problem: find a point x € fi and the smallest radius r > (if they exist) 
such that the set x + rF intersects all of the target sets. It is obvious that when F is the 
Euclidean closed unit ball of the plane, the target sets under consideration are singletons, 
and the constraint set is the whole plane, this problem reduces to the smallest enclosing 
circle problem. 

The so-called minimal time function below plays a crucial role in the study of the smallest 
intersecting ball problem. For a nonempty subset QCI, define 

Tn(x)=mf{t>0:(x + tF)nQ^$}, (1) 
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which signifies the minimal time for the point x to reach the target set Q following the 
dynamic F. General and generalized differentiation properties of the minimal time function 
(1) have been studied extensively in the literature; see, e.g. [5, 6, 8, 13, 15] and the references 
therein. When F is the closed unit ball of X, the minimal time function becomes the familiar 
distance function 

d(x; Q) = inf{||x - q\\ : q £ Q}. 

In [4-6, 8], the minimal time function (1) in which F is a closed bounded convex set that 
contains the origin as an interior point is considered. Further extensions to the case where 
the origin is not necessarily an interior point of F have been considered in [10, 13]. However, 
to the best of our knowledge, the minimal time function (1) with F being unbounded has 
not been addressed in the literature. In this paper, we are going to initiate the study of the 
minimal time function (1) for the case where F is unbounded and consider applications to 
the corresponding extended version of the smallest intersecting ball problem. 

Throughout the paper, we make the following standing assumption, unless otherwise 
specified: 

X is a normed space and F C X is a nonempty closed convex set. 

Under natural assumptions, the smallest intersecting ball problem generated by the 
dynamic F with the target sets Oj for i = 1, . . . , m and the constraint set Q can be modeled 
as the following optimization problem: 

minimize T(x) subject to x G fl, (2) 

where 

T{x) = max{7^(x) : i = 1, . . . , m}. 

The functions involved in this problem are nondifferentiable with extended-real values 
in general, and they are convex when all of the sets therein are convex. Thus it is natural 
to use nonsmooth analysis and convex analysis in particular to study the problem. 

Given a convex function <p : X — > (— oo, oo] and a point x G dom <p, the subdifferential 
in the sense of convex analysis of the function <p at x is defined by 

dip(x) = {x* G X* : (x*, x — x) < <p(x) — <p(x) for all x £ X}, 

where the set X* is the dual of X and the domain of the function tp is 

dom <p = {x G X : (f(x) < oo}. 

An element x* G dip(x) is called a subgradient of the function tp at x. It is well known that 
subdifferentials of convex functions play an important role in convex optimization in both 
theoretical and numerical aspects. For instance, the generalized Fermat rule 

ip has an absolute minimum at x if and only if G d<p(x) 

allows us to use subdifferentials as a tool to solve nonsmooth convex optimization problems. 
Under natural assumptions, the iterative sequence 

Xk + i = Xk — ctkxl, with a chosen x\ G X and x* k G dip(xk), 
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converges to a minimizer of the function (p; see, e.g., [1] and the references therein. 

Due to the involvement of the minimal time function (1) to the smallest intersecting 
ball problem, knowing subdifferential properties of this function provides important tools 
for the study of the problem in the convex setting. In fact, this has been shown by our 
recent papers [14, 16, 18] that generalized differentiation of the minimal time function (1) 
plays an important role in the study of set facility location problems including the smallest 
intersecting ball problem. It turns out that subdifferential properties of the minimal time 
function (1) have close connections to the asymptotic cone of F and the normal cone to Q 
at x G Q defined by 

N(x; Q) = {x* G X* : (x*, x - x) < for all x G Q}. 

In this paper, we review important properties of asymptotic cones for convex sets in 
Section 2. Section 3 is devoted to general properties of the minimal time function (1), while 
generalized differentiation of this function in the convex case is addressed in Section 4. 
Finally, in Section 5, we provide applications of the results obtained to numerically solving 
the optimization problem (2) and its connection to the smallest intersecting ball problem. 
We review the literature and provide self-contained proofs for most of the results for the 
convenience of the reader. 

Throughout the paper, X* denotes the dual space of X and (x*,x) = x*(x) is the dual 
pair between an element x* G X* and an element x G X. 

2 Asymptotic Cones 

In this section we review some facts about asymptotic cones that will be used in our results. 
For x G F, the asymptotic cone of F at x is defined by 

i^ooO) = {d G X : x + td G F for all t > 0} 
= {d G X : x + td G F for all t > 0}. 

Another equivalent definition of F oc (x) is given as 

, . F — x 

Foo(x) = n — — , 

which shows that F oc (x) is a cone and contains the origin, 0. In addition, F OQ (x) is closed 
and convex since F is a closed convex set and the intersection of closed (convex) sets is 
closed (convex). Thus -Foo^) is a closed convex cone. 

Proposition 2.1 For all Xi,x% G F, we have 

Foo{x\) = F 00 (x 2 ), 

that is, the asymptotic cone does not depend on x G F. Thus we will denote the asymptotic 
cone as F^. 



3 



Proof. For the proof of this statement in finite dimensions, the reader is referred [9]. 
The proof remains valid in infinite dimensions, and we include the details here for the 
convenience of the reader. Fix x\ and x 2 in F. It suffices to show that 

Foo (xi) C F OQ (x 2 ). 

Let d G F OQ (xi). Fix any t > 0, we will show that x 2 + td G F. Consider the sequence 

x n = xi+td+ ( 1 - - J (x 2 - xi), 



which can be written as 



x n = -(xi+ ntd) +1 ) x 2 . 

n \ n 



We see that x n G F for every n because d G F^xi) and F is convex. In addition, 
x n — > X2 + id. Thus x 2 +td £ F since F is closed. Therefore, d € ^00(^2)- D 
The following corollary is a direct consequence of Proposition 2.1. 

Corollary 2.1 Suppose F contains the origin. Then 

F>c = f| *F. 

i>0 

Proof. Through the use of definitions of F^ and Proposition 2.1, we have 

F oo = F oo (0) = f]^ = f]lF=f]tF 

where t = 1/i for t > 0. □ 

Proposition 2.2 The following are equivalent: 

(i) deFoo. 

(ii) There exists a sequence {t n } C [0, 00) s-uc/t i/tai i n — > and a sequence {f n } ^ F wit/i 

tnfn ^ 

(iii) There exists a sequence {t n } C [0, 00) suc/i ffoai i n — > and a sequence {fn} Q F with 

tnfn ^ d. 

Proof. Let us prove "(i)=^ (ii)"- Suppose d £ F^ and fix S G F. Then 

x + G F for all t > 0. 

In particular, 

x + nd G F for all n G IV. 
For each n, there exists f n £ F such that 

x + nd = f n or equivalently — x + d = —f n - 

n n 

Let t n = —. We easily see that t n f n — > d as n — > 00. 
n 
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The implication "(ii)=> (iii)" is obvious since the strong convergence in X implies the 
weak convergence. We finally prove "(iii)=> (i)". Assume that there exist sequences {t n } C 
[0, oo) and {f n } Q F such that t n — > and t n f n — >■ d. Fix any i£f, and let t > 0. We 
will prove that d € by showing that 

x + td e F for all t > 0. 

One has < t ■ t n < 1 when n is sufficiently large and t ■ t n — > 0. Thus 

(1 — t ■ t n )x + t ■ t n f n x + td as n —> oo. 

The elements (1 — t ■ t n )x + t ■ t n f n are in F because F is convex, and since F is weakly 
closed, x + td G F. Therefore, d € F^. □ 

Definition 2.1 T/ie set F is called sequentially compact at if for any sequence {/«} C F 
such that f n ^> 0, one has f n — ^ 0. 

This property holds obviously in finite dimensions. Moreover, if X is a Hilbert space and F 
is a nonempty closed convex subset of the Frechet normal cone to a set that is sequentially 
normally compact (see [12]) at the reference point, then F is automatically sequentially 
compact at 0. In infinite dimensions, it is not applicable to sets with nonempty interior. 

Proposition 2.3 If F is bounded, then F^ = {0}. The converse holds if we suppose further 
that X is a reflexive Banach space and F is sequentially compact at the origin with £ F 
or X is finite dimensional. 

Proof. Suppose F is bounded and d G F^. By Proposition 2.2, there exist sequences {t n } 
such that t n — > + and {f n } C F such that 

t n fn — > d as n — > oo. 

Since F is bounded, the sequence {t n f n } converges to 0, so d = 0. 

To prove the converse, we suppose by contradiction that F is unbounded and F^ = {0}, 
while F is sequentially compact at the origin with € F or X is finite dimensional. Then 
there exits a sequence {x n } C F with ||x n || — > oo. The sequence defined by 

, %n 
«n = -n TT 



is bounded. So it has a subsequence (without relabeling) such that d n —> d € F. Suppose 
first that F is sequentially compact at the origin with £ F. Since G F and -F is convex, 
d n 6 i 7 when n is sufficiently large. The sequentially compactness of F at implies d ^ 0. 
Indeed, if d = 0, then ci n — > 0, but this is a contradiction since \\d n \\ = 1 for every n. Fix 
any and let t > 0. For sufficiently large fc, one has 

( * \ * 

u k = 1 - n TT x + Ti u x k e 

V \\ x k\\J iFfcll 

since F 1 is convex. Since ■Ufc — >• a; + and F is weakly closed, x + td £ F 1 . Thus d G Foo, 
which is a contradiction. The proof is similar in the case X is finite dimensional. □ 
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3 General Properties of Minimal Time Functions with Un- 
bounded Dynamics 

In this section we study general properties of the minimal time function (1) that will be 
used in the next sections. These properties are also of independent interest. 

It is obvious that if x £ Q, then Tq(x) = 0. Moreover, Tq(x) is finite if there exists 
t > such that 

(x + tF) n Q ^ 0. 

This holds in particular when £ hit F. If no such t exists, then Tq(x) = oo. Thus the 
minimal time function (1) is an extended real-valued function in general. 

Theorem 3.1 Suppose one of the following: 

(i) X is a reflexive Banach space, Q is a bounded weakly closed set, and F contains the 
origin 0. 

(ii) Q is a compact set and F contains the origin 0. 
Then 

Tq(x) = if and only if x £ Q — F^. 

Proof. Let us prove the theorem under assumption (i). Suppose 7q(x) = 0. Then for 
every n E JN, there exists t n > such that 

(x + t n F) n Q + 0, 

and {tn} converges to 0. Then there exist u; n € Q and f n ^F for every n € IN with 

% ~\~ tnfn — ^n- 

Since Q is bounded and weakly closed in a reflexive Banach space, there exists a subsequence 
(without relabeling) such that the sequence {w n } converges weakly to lu G Q. Thus 

tnfn — >W-iaSIl-> OO. 

It follows from Proposition 2.2 that lu — x € i^, and hence x £ Q — F^. 

Let us prove the opposite implication. Fix any x £ Q — F^. Then there exist uj € Q 
and d £ -Fqo such that 

x = to — d. 

Since £ F and d £ F^, then by Corollary 2.1, n{to — x) = nd £ F for all n £ IN. Thus 

io - x £ -F for all n £ M. 

n 

1 1 
This implies (x + -F) n Q ^ for all n £ IN. Thus < Tn(x) < - for all n £ IN. 

n v n 

Therefore, Tq(x) = 0. The proof of the theorem under assumption (ii) is similar. □ 

Remark 3.1 The proof of Theorem 3.1 shows that, if Q is a nonempty set (not necessarily 
bounded) and £ F, then Tq(x) = whenever x £ Q — F^. 
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The following example is a demonstration of Theorem 3.1. 

Example 3.1 Let F = M x [—1,1] C M 2 and let Q be the disk with center at (1,0) and 
radius r = 1. Then using the definition of Fqo in Corollary 2.1 we obtain 

F 00 = Mx {0}, 

and Q - Foo = Mx [-1, 1]. Thus 7j(x) = if and only if x € -K x [-1, 1]. 

Let (p : X —> (—00,00] be an extended real- valued function and let x G dom 93. Recall 
that <p is lower semicontinuous at x if 

liminf > ^(x). 

X^tX 

The function ip is called lower semicontinuous if it is lower semicontinuous at every point 
of its domain. 



Theorem 3.2 In the same setting of Theorem 3.1, the function Tq is lower 



semicontinu- 



ous. 



Proof. Fix any x G dom Tq and fix a sequence {x n } that converges to x. We will prove 
that 

liminf Tq (x n ) > Tq{x) 

under the assumption (i) of Theorem 3.1. The inequality holds obviously if liminf Tq (x n ) = 
00. So we only need to consider the case where liminf Tq (x n ) = 7 G [0, 00). It suffices to 
show that 7 > Tq(x). Without loss of generality, we can assume that Tq \x n ) — > 7. By the 
definition of the minimal time function (1), there exists a sequence {t n } C [0, 00) such that 

Tq (x n ) <t n <TQ (x n ) + - and (x n + t n F) n Q + for all n € IN. 
For every n € IN, fix f n G F and u n G Q such that 

Xn tnfn — ^ni 



and we can assume without loss of generality that oj n ^> u) G Q. We will consider two cases: 
7 > and 7 = 0. If 7 > 0, then 

U^^eF. 

7 

Thus (x + 7F) n Q ^ 0, and hence i>Tq(x). 

Now we consider the case where 7 = 0. In this case the sequence {t n } converges to 0, 

and 

l-njn ' X. 

By Proposition 2.2, ui — x G F^. This implies x G — i 7 ^. Employing Theorem 3.1, 
one has that Tq(x) = and we also have 7 > Tq(x). We have showed that Tq is 
lower semicontinuous under the assumption (i) of Theorem 3.1. The proof for the lower 
semicontinuity under the assumption (ii) is similar. □ 
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Theorem 3.3 Suppose Q C X is a nonempty convex set. Then the minimal time function 
(1) is convex. 

Proof. Fix any x±,x 2 G dom Tq and let A G (0, 1). We will show that 

T F {Xxi + (1 _ A)X2 ) < XT F {xi) + (1 _ A )T^(X 2 ). (1) 

Let 7i = Tq(xi) for i = 1,2. Given any e > 0, there exist ii for i = 1, 2 with 

li <ti<ji + e and (xj + UF) n Q ^ 0. 
Since both F and Q are convex, one has 

[Axi + (1 - A)x 2 + (Ati + (1 - A)i 2 )F] n Q ^ 0. 

This implies 

T(f (Axi + (1 - A)x 2 ) < Ati + (1 - X)t 2 < AT^(xi) + (1 - A)T^(x 2 ) + e. 

Since e is arbitrary, (1) holds and the proof is now complete. □ 
The Minkowski function associated with F is defined by 

p F (x) = inf{i > : x G tF}, x G X. (2) 

In contrast to the familiar definition of the Minkowski function, in this definition we only 
assume that F is a nonempty closed convex subset of X. 

Proposition 3.1 The Minkowski function pp is a positively homogenous and subadditive 
extended real-valued function. Suppose further that G F . Then pf(x) = if and only if 
G Fqq . 

Proof. We first prove that pp is subadditive, that is, 

Pf(x\ + x 2 ) < pf(xi) + Pf{x2) for all x\,xi G X. (3) 

The inequality holds obviously if the right side is infinity, that means x\ £ dom pp or 
x 2 ^ dom pp. Suppose x\,xi G dom pp. Fix any e > 0. Then there exist ti,t 2 > such 
that 

Pf{x\) < t\ < pf{xi) + e, p F (x 2 ) < t 2 < pf{x 2 ) + e, 

and 

X\ G t x F, x 2 G t 2 F. 

Since F is convex, 

Xi + x 2 G hF + t 2 F C (ti + t 2 )F. 
Thus pf{x\ + x 2 ) <t\ + t 2 < pf{x\) + pf(x 2 ) + 2e. This implies (3). 
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The fact that pp is positive homogeneous follows from the following analysis. For a > 0, 
one has 

Pf{oix) = inf{i > : ax € tF} 

= inttt >0:i£ — F} 
a 

= ainf{— >0:i£ — F} = app(x). 
a a 

Now, we will prove the second statement. Let Q = {0}. We have Pf(x) = Tq F {x). Since 
€ F, by Theorem 3.1, Pf( x ) = if and only if x € Q — (—F)^ = F^. The proof is now 
complete. □ 

Recall that a function ip : X — > M is said to be £— Lipschitz if 

— p(y)\ < £\\x — y\\ for all x,y & X. 
Proposition 3.2 Suppose 6 mi -F. Define 

£ = inf{- : JB(0;r) C F,r > 0}. 
r 

TTien pf is an £— Lipschitz function. In particular, pf(x) < £\\x\\ for all x € X. 

Proof. It is clear that £ is a nonnegative real number since € int F. Let r > satisfy 
2B(0;r) C F. Then 

p F (a;) = inf{t > : x € tF} < mf{t >0:iE iiB(0;r)} = inf{t > : 11x1 1 < rt} = — . 

r 

This implies Pf( x ) < ^IMI- Since pf is subadditive, for any x,y £ X, one has 

p F (s) - p F (y) < p F (a; - y) < £\\x - y\\. 

This implies the £— Lipschitz property of pp- □ 

Using the definitions of the minimal time function (1) and the Minkowski function (2), 
we can establish a relationship between these functions in the proposition below. 

Proposition 3.3 Let QCI be a nonempty closed set. Then 

T Q F (x) = m£{p F (u - x) : u € Q}. (4) 

Proof. Let us first consider the case where Tq(x) = oo. Then the set 

{t > : (x + tF) n Q + 0} = 0. 

It follows that for every uj £ Q, the set 

{t > : w - x G iF} = 0. 

Thus /3_f(o; — x) = oo for all uj E Q, and hence the right side of (4) is also infinity. 
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Now, suppose Tq(x) < oo. Fix any t > such that 

(x + tF)nQ^0. 

Then there exist / € F and uj € Q with x + tf = uj. Thus a; — a; £ ii 7 , and hence 
Pf{u — x) <t. This implies 

inf p_p(w — ar) < t. 



It follows by the definition of the minimal time function (1) that 

'q 



inf pf(u — x) < Tq(x) < oo. 



Let 7 = inf^gQ pp{oj — x). Then for any e > 0, there exists cu £ Q with 

/9f(oj - x) < 7 + e. 

Employing the definition of the Minkowski function (2), one finds t £ [0, oo) with 

iv — x £ tF and £ < 7 + e. 

Thus wGii IF, and hence 

Tq (x) < i < 7 + e. 

Since e is arbitrary, one has 

7q(x) < 7 = inf - x). 

The proof is now complete. □ 

Employing the £— Lipschitz property of the Minkowski function (2) from Proposition 3.2 
and the representation of the minimal time function (1) in Proposition 3.3, we obtain the 
result below. 



Theorem 3.4 Let £ be defined in Proposition 3.2. Suppose £ int F . Then the minimal 
'Q 



time function Tq is £— Lipschitz. 



Proof. Let x,y £ X. Then using the subadditivity of pf from Proposition 3.1 and using 
Proposition 3.2, one has 

Pf{oj - x) < pf(u - y) + Pf(u — x)< pf{u — y)+£\\y — x\\ for all u> £ Q. 

It follows from Proposition 3.3 that 

T£{x)<T£(y)+t\\x-y\\. 

This implies the £— Lipschitz property of Tq . □ 
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4 Convex Subdifferentials of Minimal Time Functions with 
Unbounded Dynamics 

In this section, we discuss properties of the convex subdifferentials of the minimal time 
function (1) where the target set Q is convex. In this case the minimal time function (1) is 
also convex by Theorem 3.3. Consider the following extended real-valued function defined 
on X*: 

a F (x*) = sup{(x*, /) : / € F}. 

Theorem 4.1 Consider the minimal time function (1) with a nonempty convex target set 
Q. For x G Q, one has 

dT£(x) = N(x;Q)nC*, 

where 

C* = {x* G X* : a F {-x*) < 1}. 
Suppose further that G F. Then 

dT^(x) = N(x; Q-F OQ )nC* = N(x; Q)(~)C*. 

Proof. Fix any x* G 8Tq(x). One has 

(x*,x-x) <Tq(x)-Tq(x) for allxGX (5) 

Since Tq(x) = for all x G Q, one has 

(x*,x — x) < for all x G Q. 

This implies x* G N(x; Q). For any / G F and t > 0, it follows from (5) that 

(x*,(x-tf)-x) <r£(x-tf) <t. 

Note that the last inequality holds since {{x — tf) + tF) flQ^0. Thus 

<**,-/)< 1 for all / G F, 

so x* G C*. Thus dT£ (x) C N(x; Q) n C*. 
Now, we will show that 

N{x;Q)nC* CdTg(x). 
Fix any x* G N(x; Q)nC*. Then (x*, -f) < 1 for all / G F and 

(x*,x — x) < for all x G Q. 

Fix any u G dom Tq ■ For any e > 0, there exist t G [0, oo), / G i 7 , and w G Q with 

7q(«) < i < Tq{u) + e and u + tf = to. 



11 



Then 



(x* , u — x) = (x* ,u — oj) + (x* ,U) — x) 

< (x*, -tf) < t < Tq(u) + e = T£{u) - T£(x) + s. 

This implies (5) since e is arbitrary. 
Let us finally prove that 

dT$ (x) = N(x; Q - Foo) n C* = N(x; Q) n C* , 

under the additional assumption. Obviously, N(x;Q — F^) n C* C N(x;Q) n C* since 
Q C Q - Fqo. Fix x* G iV(x; Q) (~l C* and any x G Q - F^. One has x = w- (iforu;GQ 
and d 6 Fx,. Since G F, one has td & F for all £ > 0. This implies (sc*, —id) < 1 for all 
i > 0. Thus (x*, -d) < 0. Then 

(x* , x — x) = (x* ,u — d — x) = (x* ,u — x) + (x* ,—d) < 0. 

Therefore, x* £ N(x; Q — F^) n C*. The proof is now complete. □ 

Let us now consider the case where x G Q — F^. In fact, similar results like those in 
Theorem 4.1 hold true in this case. 

Theorem 4.2 Consider the minimal time function (1) with a nonempty convex target set 
Q. Suppose xGQ - Fqo and G F. Then 

dT£(x) = N(x;Q-F OD )nC*. 

Proof. Fix any x* G dTq(x). Using (5) and the fact that Tq(x) = for all x G Q — F^ 
from Remark 3.1, one has that x* G N(x; Q — F^). Since x G Q — F^, one has x = oi) — d, 
where u> G Q and d G F^. Fix any / G F. For any i > 0, let xt = {x + d) — tf = ui — tf. 
Then (x t + tF) n Q / 0. Thus 

(x*,x t - x) < ^(xt) < t, 

or equivalently, 

(x*,d-tf) < t. 

Thus 

(x*,~- f) < 1 for all i > 0. 

This implies (x*, — /) < 1 by letting t — > oo, and hence x* G C* . Thus we have shown that 
dT£(x) CiV(x;Q-F 0O )nC Y *. 

Now fix any x* G N(x; Q — F^) PI C*. Fix any u G dom Tq ■ Then there exist i G [0, oo), 
w€Q and / G F such that 

7g («) < t < Tq(u) + e and u + tf = u. 

Since Q C Q — F^, , one has 

(x*,u - x) = (x*, u - x) + t(x*, -/) < t < T£{u) + e = T£(u) - T£(x) + e. 
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This again implies x* G dTq(x), and thus N(x; Q — Fqo) n C* C &Tq(x). The proof is now 
complete. □ 

In what follows we are going to discuss the subdifferential in the sense of convex analysis 
for the minimal time function (1) when the reference point does not belong to Q — F^,. Our 
development is inspired by the early work on sub differentials of distance functions; see 
[2, 3, 11]. For any r > 0, define 

Q r := {x £ X : Tq (x) < r}. 

Lemma 4.1 Let Q be a nonempty set and let r > 0. Suppose x ^ Q r with Tq(x) < oo. 
Then 

7j(*)=7£(s)+r. 



'Qr 
'Qr 

(x + tF) n Q r + 0. 



Proof. Let r > 0. Since Q C Q r , one has T~X (x) < Tq(x) < oo. Fix any f > such that 



Then there exist fi £ F and it € Q r such that x + tfi = u. Since u G Q r , then Tq(u) < r. 
Then for any e > 0, there exists s > such that s < r + e. Consequently, there exist wGQ 
and /2 G i 7 such that n + s/2 = w. Thus 

w = u + s/ 2 = (x + tfx) + s/ 2 G x+ (t + s)F 

since F is convex. This implies 

Tq[x) < t + s < t + r + e. 
Since e > is arbitrary, Tq{x) <t + r, and hence 



T$ (x) < rl + r. 

Let us prove the opposite inequality. Let 7 = Tq(x). Then r < 7 since x ^ Q r . For 
any e > 0, there exist t G [0, 00), / G F, and w G Q with 

7 < i < 7 + £ and x + tf = u. 

One has 

w = x + £/ = x + (t - r)f + rf G x + (f - r)f + rF. 

Thus 

7^ (x + (i - r)f) < r. 

So x + (t—r)f G Q r ; in addition x+(t — r)f G x + (t — r)F. Hence Tq r {x) < t — r < 7 — r + e. 
Since e > is arbitrary, one has 

r + Tl <7 = T c f(x). 

The proof is now complete. □ 
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Lemma 4.2 Let Q be a nonempty set and let x G dom Tq ■ For t > and any f G F, one 
has 

Tf(x-tf)<T£(x)+t. 

Proof. For any e > 0, there exists s > such that 

Tq(x) < s < Tq(x) + e and (x + sF) DQ ^ 0. 

Then (x -tf + tF + sF) C[Q ^ 0, and hence (a; - tf + (t + s)F) n Q / 0. Thus 

7^(x - tf) < t + s < Tq(x) + t + e. 

The conclusion follows by letting e — > 0. □ 

Theorem 4.3 Xei X &e reflexive Banach space and let £ F. Consider the minimal time 
function (1) with a nonempty closed bounded convex target set Q. For x ^ Q — F^, one has 

dT£(x) = N(x;Q r )nS*, 

where 

S* = {x* G X* : a F (-x*) = 1} and r = Tq (x) > 0. 

Proof. Fix any x* G dTq(x). Following the proof of Theorem 4.1, see that o\f(— x*) < 1 
and x* G N(x;Q r ). Let us show that o~f{— x*) = 1. One has 

{x*,x -x) <Tq(x)-Tq (x) for all x G X. (6) 

Fix any e G (0, r). There exist t G M, f G F and wGQ such that 

r < t < r + e 2 and oj = x + tf. 

We can write w = x + e/ + (t — e)f. So Tq (x + ef) < t — e. Applying (6) for x = x + ef, 
one has 

(x*,sf) < T£(x + ef) - T£ (x) < t - e - r < e 2 - e. 

This implies 

l-e<(-x*,f) <a F (-x*). 

So a F (-x*) > 1 by letting e -»• 0. Therefore, x* G S*, and thus OTq (x) C N(x; Q r ) n 5 1 *. 

Fix any x* G A^(x; Q r ) such that a F {— x*) = 1. We are going to show that (6) is satisfied. 
By Theorem 4.1, x* G dTq r (x). Thus 

(x*, x — x) < Tq t {x) for all x G X. 

Fix any x G X. In the case t := Tq(x) > r, one has Tq{x) — r = Tq t {x) by Lemma 4.1. 
Thus (6) holds. Suppose t < r. For any e > 0, choose f £ F such that (x*, — /} > 1 — e. 
By Lemma 4.2, Tq {x — (r — t)/) < r. So x — (r — t)/ G Q r . Since x* G iV(x; one has 

(x*,x-(r-t)f -x) < 0. 

This implies 

<x*, x - x) < (x*, f)(r - 1) < (1 - e)(t - r) = (1 - e)(T Q F (x) - 7f (x)). 
Since e is arbitrary, (6) holds. Therefore, x* G dTn(x). □ 
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Theorem 4.4 Let X be reflexive Banach space and let £ F. Consider the minimal time 
function (1) with a nonempty closed bounded convex target set Q. Suppose x ^ Q — and 
€ F. Then 

9Tq(x) = [-dp F (Co - x)] n N(Q; Q) 
for any Co 6 IL F (x;Q), where 

U F (x; Q) = {u} eQ :Tq (x) = p F (oo - x)}. 
Proof. Fix any x* € dTq \x) and Co € H F (x;Q). Then Tq(x) = p F (Co — x) and 

(x*,x-x) < Tq(x) -Tq{x) for all x G X (7) 

For any oo £ Q, one has 

(x* , oj — a)) = (x* , (w — Co + x) — x) 

<T£(u:-Cd + x)-T£(x) 

< p F (uj — {oo — Co + x)) — p F (Co — x) = 0. 

This implies x* G iV(a); Q). Let u = Co — x. For any i € (0, 1) and u 6 X, applying (7) with 
x = x — — S), one has 

(x*, - u)) < T Q F (x - t(u - u)) - T Q F (x) 

< p F (Cd — (x — — u))) — p F {Cd — x) 

= p F (u + t(u - u)) - p F (u) 
= p F (tu + (1 - t)u) - p F (u) 

< tp F (u) + (1 - t)p F (u) - p F (u) 
= t{p F {u) - p F {u)). 

It follows that 

{—x*,u — u)< p F {u) — p F (u) for all u € X. 

Thus —x* € dp F (u), and hence x* € —dp F {Co — x). The inclusion "c" has been proved. 
Let us prove the opposite inclusion. We have that 

p F (x) = mf{t >0:x£tF} = inf{t > : (-x + tF) D O + 0} = Tq (~x), 

where O = {0}. For any Co € Q one has, x — Co ^ O — F^ = —Foo since x ^ Q — F^. From 
—x* £ dp F (Co — x), one has x* £ S* by Theorem 4.3. Thus we only need to show that 

[-dp F {oo - x)} n N(oo; Q) C N(x; Q r ), 

and the conclusion will follow from Theorem 4.3. 

To proceed, pick any x € Q r and for an arbitrary small e > 0. Then there exist t < r + e, 
f e F, and weQ with uo = x + tf. Thus (— x*, /) < <7p(-x*) = 1 and 

(x*, x — x) = (x*, u; — tf — x) 

= t(—x*, f) + (x* ,oo — Co) + (x* ,oo — x) 

< t + (x* ,oo — Co) + {x* ,00 — x) 

< Tq(x) + e + (x*,o; - + (x*,o) - x). 
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We have (x*,u) — ui) < due to x* G N(p; Q) and 

(x*,6j — x) = (—x*, — (ui — x)) < pp(0) — pf(ui — x) = —Tq(x) 

by —x* G 0pf(O — x). It follows therefore that {x*,x — x) < e for all x G Q r , and hence 
x* G N(x; Q r ) because e > was chosen arbitrary small. □ 

Let p> : X — > (— oo, oo] be a convex function and let x G dom p. The singular subdiffer- 
ential of p> at x is defined by 

d°°ip{x) = {x* G X* : (x*,0) G N((x,ip(x));epiip)}. 

Lemma 4.3 Let p : X — > ( — 00,00] &e a convex function and let x G dom 99. IFe always 
have 

d°°p(x) = N(x; dom tp). 
Proof. Fix any x* £ d°°<p(x) and x G dom 99. Then (x,</?(x)) G epi 93. Thus 

{x* ,x — x) = (x*,x — x) + 0(<^(x) — V'(S)) < 

because (x*,0) G N((x, p(x)); epi </?). Conversely, suppose x* G iV(x;dom<£). Fix any 
(x, A) G epi /. Then p(x) < A, so x G dom 99. Thus 

(x*, x — x) + 0(A — = (x*, x — x) < 0. 

This implies (x*,0) G N((x, p(x)); epi p) or equivalently x* G d°°ip(x). □ 
Define the following subset of X* 

F* = {x* G X* : {x* , /} > for all / G F}. 

We are going to obtain below explicit representations of singular subdifferentials of the 
minimal time function (1). 

Proposition 4.1 Consider the minimal time function (1) with a nonempty convex target 
set Q. The following hold: 

(i) If x G Q, then 

d c °T£(x) = N(x;Q)nF;. 

(ii) // x G Q — Foq and G F, then 

d°°T£(x)=N(x;Q-F oc )nF*. 

(iii) If x $l Q — F^, X is a reflexive Banach space, and Q is closed and bounded, G F, 
and Tq{x) < 00, then 

d°°T£(x) = N(x; Q r ) n F*, r = T£(x) > 0. 



16 



Proof, (i) Fix any x* £ 9°°7q \x) = N(x; dom Tq). Employing Lemma 4.3 and the fact 
that Q C dom Tq, one has d^Tqix) = N(x;dom Tq) C N(x;Q). We are going to show 
that (x*, /) > for all / £ F. For any / € F, one has x — f £ dom 7^ since Tq(x — f) < 1. 
Thus 

(a?*,(z -/) -S) < 0. 

This implies (x* , f) > 0. Let us prove the converse. Fix any x* £ N(x;Q) n F+ and 
rc G dom T^f. Then there exists t 6 [0, oo) such that (x + iF) PI Q 7^ 0. Fix further f £ F 
and oj £ Q with x + tf = w. Then 

(a;*, x — x) = (x*, u) — tf — x) 

= (x*,to-x) -t(x*,f) < 0, 

Thus we get x* £ N(x; dom 7j) = d 00 ^ 

(ii) Since Tq(x) = for all x £ Q — F^, one has Q — F^ C dom Tq ■ Applying Lemma 
4.3 again, one has 9°°7q \x) = N(x; dom Tq) C iV(x; Q — -^oo)- Let Co £ Q and d 6 Foo be 
satisfy x = u) — d. Since 6 F, one has d £ F. For any f £ F, 

x - f = oj - (d + f) C Q - 2(l/2d + 1/2/) CO-2F. 

This implies Tq(x — f) < 2, so x — / € dom 7^. This again implies x* € -F+ as in the 
proof of (i). The proof of the oppositive inclusion follows by the same proof of (i) using the 
observation that Q C Q — F^. The proof for (iii) is similar. □ 

5 Applications to the Smallest Intersecting Ball Problem 

In this section, we consider the optimization problem (2), where the target sets fij for 
i = 1, . . . ,m and the constraint set f2 are nonempty closed bounded convex subsets of X. 
The smallest intersecting ball problem (SIB) under consideration asks for a point x £ 
and the smallest t > (if they exist) such that 

(x + tF) n fii / for alH = 1, . . . , m. 

Definition 5.1 We say that x £ is a solution of the smallest intersecting ball problem 
(SIB) with the smallest radius r if r £ [0, 00), 

(x + rF) n / for all i = 1, . . . , to, 

and whenever there exist x £ O and r' £ [0, 00) with 

(x + r'F) n Qi / /or all i = 1, . . . ,m, 

then r < r' . 

Lemma 5.1 Let X be a reflexive Banach space. Suppose Q C X is a nonempty closed 
bounded convex set and 6 F. IfT£(x)<r and r > 0, then (x + rF)nQ ^ 0. 
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Proof. Let 7 = Tq(x). Suppose first that 7 > 0. Then there exists a sequence {tk} that 
converges to 7 such that x + tkfk = w*., where /& G F and w fc € Q. Since Q is closed convex 
and bounded, we can assume that uj^ ^ oj £ Q. Then . Since F is closed and 

convex, it is weakly closed, and hence G F. Thus (x + 7F) n Q 7^ 0. It follows that 

7 

(x + rF) (~l Q 7^ since 7 < r and G F. If 7 = 0, then 1 £ Q - F^, by Theorem 3.1. Let 
x = oj — d, where uj £ Q and d G F^. By Corollary 2.1, for r > 

w — x d „ 
= - G F. 

r r 

This again implies (x + rF) n Q 7^ 0. The proof is now complete. □ 

Theorem 5.1 Let X be a reflexive Banach space and let G int F. Suppose £1 and for 
i = 1, . . . , m are nonempty closed bounded convex sets and 

on[n^ 1 (n i -F oo )] = 0. (8) 

TTien x £ £1 is an optimal solution of the optimization problem (2) with r = T(x) if and only 
if x is an optimal solution of the smallest intersecting ball problem (SIB) with the smallest 
radius r. 

Proof. Suppose x G is an optimal solution of (2) and r = T(x). If r = 0, then T£.(x) = 
for all i = 1, . . . , m. This implies aiGfij- for all z = 1, . . . , m, which is a contradiction 
to (8). Thus r > 0. From Lemma 5.1, (x + rF) n fij 7^ for all i = 1, . . . , m. Let us fix an 
x G f2 and r' > such that (x + r'F) H f2j 7^ for alH = 1, . . . , m. We will show that r' > r. 
Suppose r' < r. Since 7^(x) < r 1 for all i = 1, . . . , m, one has T{x) < r' < r = T(x), 
which is a contradiction. 

Let us prove the converse. Suppose that x is an optimal solution of the smallest inter- 
secting ball problem (SIB) with the smallest radius r. From Definition 5.1, one has x £ 0, 
and (x + rF) fl Oj 7^ for all i = 1, . . . , m. Thus 

T{x) = max{7^(:r) : i = 1, . . . , m} < r. (9) 

By (8), r > 0. Fix any x G £1 and define r' = T(x). Since G int F and by (8), one 
has r' G (0, 00). By the definition of the minimal time function (1), 7~q.(x) < r' for all 
i = 1, . . . , m. It follows again from Lemma 5.1 that 

(x + r'F) n Sli ^ for alH = 1, . . . , m. 

Thus r < r' = T{x) by Definition 5.1. In particular, taking x = x, one has r < l~(x), 
and hence r = l~(x) by (9). So T(x) < T{x). Therefore, rr in an optimal solution of the 
optimization problem (2). □ 

For each u G X, define 

/(u) = {»€{!,..., m}:T(«)=7£(u)}. 
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Theorem 5.2 Let X = M n with the Euclidean norm and let G int F. Consider the 
optimization problem (2) in which for i = 1, . . . ,m and are nonempty closed bounded 
convex sets. Fix x\ G and define the sequences of iterates by 



where 



xl G < 



^fc+i = V(x k - a k x* k ; O), Vfc G IV, 

{0} i/ x k G fij 

[-9pp(a; fc - x k )\ n iV(w fc ; Oj) i/ x k £ 



for some i G I(x k ), and V denotes the Euclidean projection operator. 

Let V k = mm{F(xj) : j = l,...,k}. If Y^h=i a i = 00 an d a l < °°> then the 

sequence {x k } converges to an optimal solution of problem (2) and V k converges to the 
optimal value V of the problem. Moreover, 



Vk _y<^- Xl]S)2 + L2E ^ a * 



2£f=i 



a, 



where < L < oo is a Lipschitz constant of the function T, and S is the solution set of 
problem (2). 

Proof. First observe that 5^0 since T is a continuous function and the constraint Q 
is compact. We have the following well-known formula, see e.g. [9, Corollary 4.3.2], for 
computing the subdifferential of the function T at a point u: 

dT(u)=co{dTi i (u):ieI(u)}. 

Thus for any i G L(u), one has dT^{u) C dT{u). The choice of x* k in the statement of 
the theorem allows us to find an element of dT^{x k ) by Theorem 4.1, Theorem 4.2, and 
Theorem 4.4. Then x k G dT{x k ). The rest of the proof comes from the standard projective 
subgradient method; see [1]. □ 

In the proposition below, we provide explicit formulas for computing the asymptotic 
cone and the Minkowski function generated by an unbounded dynamic. 

Proposition 5.1 Let x* G X* and 6j > for i = 1, . . . , q. Define 

F = {x G X : {x*,x) < bi for all i = l,...,q}. (10) 

Then 

Foo = {x G X : (x*,x) < for all i = 1, . . . , q} 

and 

(x* , u) 

Pf(u) = max {0, max{ — j- — : i = 1, . . . , g}|. 
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Proof. Let x* G X* . By the definition of F we have that F is a closed convex set. Since 
£ -F, by Corollary 2.1, one has that ci G if and only if 

td 6 F for all t > 0, 

or 

(x*,d) < -j for all i = 1, . . . , q and for all t > 0. 

This is equivalent to 

(a;*, d) < for all i = 1, . . . , q. 

The formula for F^ has been verified. 

Since bi > for all i = 1, . . . , q, one has G int -F, and 

p F (u) = M{t >0:uetF} 

= inf{t > : (x*,u) < tbi for all i = 1, . . . , q} 

= mf{t > : ^h^t < t for alH = 1, . . . , q} 
= max { 0, max{ — - : i = 1, . . . , g}} . 

The proof is now complete. □ 

The following corollary provides a formula for computing the subdifferential of Pf{ u ) 
when u £ F^. The formula allows us to determine a subgradient x£ in Theorem 5.2 when 
F is the polyhedron dynamic (10) in the case Xk ^ &i — F^. 

Corollary 5.1 Consider the set F defined in Proposition 5.1. For u G X, define 

J(u) = {i = l,...,p:^- = p F (u)}. 

For any u £ Fx,, one has J(u) ^ and 

dp F {u) = co :ie J(u)}. (11) 
bi 

Proof. Since u £ F^, by Proposition 5.1, there exists j = l,...,q such that (x*,u) > 0. 
This implies Pf(u) > and J(u) ^ 0. Using the formula for computing pp in Proposition 
5.1 and the well-known formula for computing subdifferential of "max" function, see e.g. 
[9, Corollary 4.3.2] (which also holds in this setting), one obtains (11). □ 

The following example shows an application of the subgradient method. 

Example 5.1 Let F = {(x, y) G M 2 : x < 1, y < 1} where x\ = (1,0) and x\ = (0,1). 
Then F^, = {(x,y) : x < 0, y < 0}. Define the target sets fij by two-tuples, that is, 
fli = {(di, bi)} for i = 1, . . . , m and let Q, = S3(c; r). Then p F ((ui, 112)) = max{ui, U2, 0}, 
and 7jf ((1*1,1*2)) = PF{(o<i,bi) — (1*1,1*2)) = max{0, aj — i*i,£>j — U2}. Suppose that we 
have found Xk = (xik,X2k) from the subgradient algorithm. Then we compute T^(xk) = 
max{0, di — xik,h — X2k} and T(xk) = max{7^(xfc) : i = l,...,m}. Fix i G I(xk). If 
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x ik > a i and X2k > bi, then choose = 0. Otherwise, x^ ^ — -Foo- If «i — ^lfc > — ^2fci 
then choose = (—1,0). In the case ai — xi k < bi — x 2k , choose x* k = (0,-1). The 
Euclidean projection to Vt can be easily obtained. Thus we can determine x k +i, and the 
algorithm is now explicit. 
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Figure 1: A smallest intersecting ball problem to squares with an unbounded dynamic. 

Example 5.2 Let F be the same set as in the previous example. Let be a square of 
right position [con — ri,uu + rj] x [a;2i — ^i-,^2i + fi\ with center Cj = (uu,uj2i) and short 
radius rj for i = 1, . . . , m. Define the constraint ft to be a nonempty closed bounded convex 
set. In this case, 

&i - Foe = {(xi,x 2 ) : x\ > oju - Ti,x 2 > uJ2i - n} 

and the projection Q can be represented by Cj and r^. 

The sequence {x* k } in Theorem 5.2 can be determined as follows 

(-1,0), if xi k < LJu - Ti and x 2k - u 2i > x lk - wi„ 

(0, -1), if x 2k < u 2 i - Ti and x 2k - uj 2i < x lk - u>u, 

(0,0), otherwise 

where i € I(xk)- The sequence {Vk} can be determined based on the fact that i 6 I(xk), 
so we have 



T ( (A(x lk ,x 2k )) = < 



(oju - n) - x^, 

(u 2l -Ti)- X 2k , 

0. 



if xik < 0Ju - Ti and x 2 k - u 2l > x\ k - u>u, 
if x 2k < lo 2i - n and x 2k - lq 2i < x± k - uiu, 
otherwise. 



Consider the target sets fij, i = 1, 2, . . . , 8, to be the right squares with centers (—30, 0), 
(-26,-34), (-20,-15), (-13,8), (-9,-25), (2,-1), (8,-30) and (8,5) with radii n = 2 
for all i = 1, 2, . . . , 8, and the constraint set is the disk with center (—7, —4) and radius 4. 
Choose a k = 1/k that satisfies Theorem 5.2 and the starting point (—7,-4). An optimal 
solution and the optimal value are x « (-3.10208, -3.10208) and V « 9.10208; see Figure 1. 
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6 Concluding Remarks 



This paper is a part of our project involving set facility location problems. The main 
idea is to consider a much broader situation where singletons in the classical models of 
facility location problems are replaced by sets. After developing new results for general and 
generalized differentiation properties of minimal time functions with unbounded dynamics, 
we are able to provide new applications to the smallest intersecting ball problem with convex 
target sets. 

Our next goal is to study the smallest intersecting ball problem with unbounded dy- 
namics and nonconvex target sets. To achieve this goal, we need to develop nonconvex 
generalized differentiation theory for minimal time functions generated by unbounded dy- 
namics. We foresee the potential of success of this future work. 
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